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Abstract 

We investigate density fluctuations in a scenario with gravitino dark matter 
in the framework of modulated reheating, which is known to generate large non- 
Gaussianity. We show that gravitino dark matter is disfavored in this framework. 
We also briefly discuss the case with the curvaton mechanism and some other possible 
dark matter scenarios. 



Recent progress in high-precision observational cosmology with cosmic microwave back- 
ground (CMB) radiation has motivated us to study the tiny nonlinear effects in the very 
early stage of the Universe when the amplitude of fluctuations was much smaller than 
unity. Among various observables to quantify non-Gaussianity of primordial fluctuations, 
the so-called non-linearity parameter /nl, which is a stunted parameter to quantify the 
degrees of freedom of actual non-Gaussianity, has been widely discussed recently because 
observations of CMB is becoming to be able to measure its possible deviation from zero. 
From the recent result from WMAP5, the so-called local- type non-linearity parameter /nl 
is constrained as —9 < /nl < HI [1] and —4 < /nl < 80 [2] at 95 % confidence level 
(CL). Although purely Gaussian fluctuation with /nl = is perfectly consistent with the 
data, its observed central value is so much deviated from zero that future higher-precision 
observations may well confirm this local-type non-Gaussianity. Then the simplest models 
of single-field inflation [3] would be ruled out as a generation mechanism of primordial 
density /curvature fluctuations [4]. 

We can think of two mechanisms as an alternative to the inflaton to generate primor- 
dial fluctuations, modulated reheating [5] and the curvaton [6,7], which are motivated in 
particle physics because there are many light fields in generic supersymmetric theories. 
Neither of them is necessary as long as the inflaton can be responsible for the generation 
of primordial fluctuations. But if the presence of large non-Gaussianity was established, 



they would be a target of serious studyr H In fact, non-Gaussianity in these alternative 
scenarios has been studied already in [9,10] for modulated reheating and in [11,12] for the 
curvaton, and it has been shown that in both scenarios large local-type non-Gaussianity 
is possible depending on the values of model parameters and/or initial conditions. 

The purpose of this Letter is to argue that confirmation of non-Gaussianity of local- 
type would not only rule out the simplest single-field slow-roll inflation models but also 
constrain the identity of the dark matter. That is, we show that the gravitino dark matter 
would confront difficulty in the two representative generation mechanisms of fluctuations 
mentioned above. Before proceeding to this main topic, for completeness, let us mention 
the relevance of the other type of non-Gaussianity, namely the equilateral-type, to the 
nature of the dark matter. It is generated due to nonlinear interactions during inflation. 
Typical models generating such non-Gaussianity are those making use of scalar fields with 
non-canonical kinetic terms [13]. They do not give us any information on the nature 
of the dark matter, because the inflaton is responsible for the fluctuation in this case. 
Observationally, the non-linearity parameter of the equilateral- type, /nx" 1 , has been much 
less constrained than the local-type, and as long as the former is concerned, there is no 
hint of non-Gaussianity so far. We therefore concentrate on the local-type non-Gaussianity 
which we denote simply by /nl as above. 

In models with a supersymmetric extension of the standard model, the gravitino, the 
spin-3/2 superpartner of the graviton, arises when the global supersymmetry is promoted 
to a local symmetry. If it is the lightest supersymmetric particle (LSP) and stable be- 



Other phenomenological generation mechanisms of large local-type non-Gaussianity include Refs. [8]. 
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cause of the .R-parity conservation, it can be a good candidate of dark matter (DM) of the 
Universe. Gravitinos can be produced in scattering processes of particles in the thermal 
bath and its thermal relic density depends on the efficiency of the production, which is 
sensitive to the reheating temperature Tr after inflation. Its relic abundance is given as 
Y3/2 = n 3 / 2 / s Tr where Y" 3 / 2 is the yield from the thermal production, n.3/2 and s are 
the number density of gravitinos and the entropy of the universe [14]. Recently, it was 
pointed out that gravitinos are also produced non-thermally from the decay of heavy scalar 
fields like an inflaton or moduli if they acquire non-vanishing expectation values [15]. The 
abundance of gravitinos produced non-thermally also depends on the reheating temper- 
ature, in fact, is proportional to the inverse of the reheating temperature. Then, when 
one considers the modulated reheating scenario as a generation mechanism of primordial 
fluctuations, the dependence on T R would result in the generation of isocurvature fluctua- 
tions in gravitino DM since the modulated reheating scenario indicates that 5Tr/Tr 7^ 0. 
Such isocurvature fluctuations can be strongly constrained from cosmological observations 
such as CMB. From WMAP5, the ratio of isocurvature fluctuations to curvature fluctua- 
tions, which is denoted by a, is constrained as a < 0.072 (axion type; uncorrelated) and 
a < 0.0041 (curvaton type; anti-correlated) for WMAP+BAO+SN at the 95% CL [1] (see 
also Ref. [16] for the constraint on isocurvature fluctuations from the non-Gaussianity). 

One may classify the origin of gravitino DM to the following two cases in the context 
of the modulated or the curvaton scenario. The first one is that gravitino can be produced 
at the decay of the inflaton through the rescattering process of thermal plasma (thermal 
production) or directly by the inflaton decay (non-thermal production). In this case, such 
gravitino DM results in the generation of large isocurvature fluctuations not only in the 
modulated scenario but also in the curvaton scenario [11]. The other one is the case 
that the curvaton mechanism is operative and gravitino DM is mostly produced from the 
decay of the curvaton. In this case, isocurvature fluctuations are suppressed if and only 
if the curvaton practically dominates the energy density of the Universe just before its 
decay [11]. However, such a case cannot generate the large non-Gaussianity. Thus when 
large non-Gaussianity is confirmed in the future, a scenario with gravitino DM may not 
be a viable candidate as the main constituent of the DM in the Universe. 

In this Letter, we discuss in detail the viability of the gravitino DM scenario when 
non-Gaussianity is large, particularly focusing on the modulated reheating scenario [5] 
since the issues of isocurvature fluctuations in this scenario is relatively unexplored. As 
shown below, the gravitino DM acquires large isocurvature fluctuations in this scenario, 
which is too large to be consistent with current observations of CMB. 

Now we briefly review the non-linearity parameter which characterizes non-Gaussianity 
of the fluctuations. For this purpose, it is fashionable to write the comoving curvature 
perturbation, up to the second order to discuss the bispectrum as 



where is the Gaussian part of curvature perturbation and /nl characterizes the size of 
non-Gaussianity observable in the bispectrum. With this parameterization, the definition 



C = C(l) + e/nlC(1) , 
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of the bispectrum is given by 

(CfcCfeCs.) = (^fBdku h, h)5(h + k 2 + k 3 ). (2) 
The leading contribution of the bispectrum can be written as 

Bc(k 1 ,k 2 ,k 3 ) = p^Pdk^ih) + P ( (k 2 )P c (k 3 ) +P c (fc 3 )P c (fc 1 )). (3) 
Here P^ is the power spectrum which is defined as 

(CfeCft,> = (27r) 3 n(W*i + fe)- ( 4 ) 

In order to evaluate the curvature perturbations (, we adopt the 5N formalism [17]. 
The curvature perturbation on sufficiently large scales ( at time tf is equal to the perturba- 
tion in time integral of the local expansion from an initial flat hypersurface (t = tj) to the 
final uniform energy density hypersurface. Since the local expansion on sufficiently large 
scales can be approximated by the expansion of the unperturbed Friedmann universe, the 
curvature perturbation ( is evaluated as 

((tf,x) = N(ti,tf,x) — (spatial average), (5) 

where the number of e-folds N(U,tf, x) is defined by the time integral of the local Hubble 
parameter as, 

ftf 

N(U,t f ,x) = / H(t,x)dt. (6) 
Jti 

In the SN formalism, taking ti = which is a certain time well after the relevant scale 
crossed the horizon scale, the curvature fluctuations from fluctuations of scalar fields are 
given as, 

t = N a 5 V :+ l -N ab 5^ a M + --- ) (7) 

up to the second order, where §(p® is the long-wave frozen part of the fluctuation of the 
scalar field (p a on the initial flat hypersurface at t = t* and N a and N a b are given by 

ON Ar d 2 N 



dip a d(p a dLf b 

Here the summation is implied for the repeated indices. 
The non-linearity parameter is then given by 



6 N ab N a N b 

(10) 



r/jVL — ,° r 2 ' (9) 



where the indices are lowered and raised by using the Kronecker's delta S ab . 

Using 5N formalism, we can evaluate the curvature perturbations in the modulated 
reheating scenario. In this scenario, fluctuations, 8a, of some modulus, a, renders fluctua- 
tions in the decay rate of the inflaton. To evaluate £ from the fluctuations of the modulus, 
we need to follow the background evolution from the epoch after the inflaton begins to 
oscillate around the minimum of the potential whose shape is assumed to be quadratic 
below. Then the energy density of the oscillating inflaton field behaves like matter, thus 
the background equations after inflation are governed by the following equations: 

r 

3 P<t> = —EP4>i (11) 



dN rv H' 

3* + = (12) 

H 2 = \{pt + p r ), (13) 

where Y is the cr-dependent decay rate of the inflaton and p^ and p r are energy densities 
of inflaton and radiation, respectively. We assume the energy density of the modulus is 
negligible throughout. The reduced Planck scale is set to be unity here. 

By solving the above equations from the end of inflation to the completion of reheating 
with the initial condition p<p(0) = po = 3iJp an d Pr(0) = 0, we can obtain a relation 
between N and T. The number of e-folds, N, elapsed while the Hubble parameter drops 
from H to Hf (Hf <C T) can be written formally as [10] 

N=ho g ^ + Q (^-), (14) 



2 Hf \ Hp 



where the function Q is defined by 
exp 



A N 'P4>( N ') 



AQ(T/H ) / dN' jjT^e™ ^-L. (15) 
J J H(N') p 

On dimensional grounds Q depends only on the combination of x = T/H . Let us consider 
the typical case that the final stage of reheating is governed by a perturbative decay of the 
inflaton after a certain period of field-oscillation dominated stage. Then x is much smaller 
than unity and we find the following analytic form of Q(x) [10] 



#2 



Q(x) = —]ogx + 0{x), x < 1. (16) 
6 

In the modulated reheating scenario, the curvature perturbation ( is evaluated as [10] 

C(N F ) = QM* + \q«MI (17) 



# 2 The behavior of the inflaton after inflation depends on the potential shape and the decay rate. 
Here we consider the case that the inflaton behaves like matter assuming the effective potential around 
the minimum is quadratic and oscillates enough before its decay. For other cases, some coefficients are 
changed but the essential results remain intact. 
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up to the second order with 



o 



Q aa = Q>( x )—^ + Q«( x )\^lj , (19) 

where Q'{x) = dQ(x)/dx and T'(a) = dT(a)/da. 

From the above equations, we can obtain the following expression for the non-linearity 
parameter 



- f NL ~ 6 1 

5 V (P(a)/r(a)) 



/nl ^6(1- ):,;7^ 2 ) (20) 



in the case i C 1. 

Hence, if we adopt proper form of the decay rate, we can realize a large non-Gaussianity 
in the modulated reheating scenario. As an example, let us consider the case that the decay 
rate depends on a as 



Via) 



1+ 4+ b (f 



(21) 



where A and B are some coefficients, M is a some energy scale and we have assumed 
a/M -C 1 and truncated at the second order in a/M. In this model, we find an simple 
expression for /nl as 

6 h L ^6(l- 2 4). (22) 



Thus, if we take |S|/A 2 3> 1 and B < 0, we can obtain a large and positive value of /nl- 
Let us move on to the isocurvature fluctuations from gravitino DM in the modulated 
reheating scenario. As mentioned before, there are several distinct ways to produce grav- 
itinos. One is the scattering of particles in the thermal plasma and their relic abundance 
is evaluated as [14] 

n/2 ^, 10 -. ?9f ( 1 + ^_)(_I|_) (23) 

in terms of the yield, where 77,3/2 is the number density of gravitino and s is the entropy 
density. In the right hand side, Tr is the reheating temperature, m 3 / 2 is the gravitino 



mass, rriGi is the gaugino masses for i-th generation, and gi is the gauge coupling^ 3 



# 3 In fact, the gauge coupling and the gaugino masses mildly depend on the reheating temperature Tr. 
However we neglect the weak dependences in the following discussions. 
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Notice that the relic abundance is proportional to the reheating temperature Tr. In the 
modulated reheating scenario, the reheating temperature fluctuates in space, which implies 
that the number density of gravitinos per entropy also fluctuates. Then, if the gravitinos 
constitute DM of the Universe, isocurvature fluctuations are generated. In order to discuss 
it more quantitatively, we relate the fluctuations of the number density of gravitinos to 
the curvature fluctuations as 

_ 5(n 3/2 /s) 

03/2 = 1 — • (24) 

n 3 /2/s 

In the modulated reheating scenario, 5Tr/Tr 7^ 0, which implies that a significant isocur- 
vature fluctuation is generated as 

^3/2 = (25) 

Since T R is proportional to T 1 / 2 , we find 8T R /T R = <5r/(2r) = Sx/(2x). Then using the 
expression for the curvature perturbation ()17j) . S'3/2 can be related to £ as 



STr 

T r 2xQ'{x) 



S3/2 = ^ = WZT^C - -3C, (26) 



when x < 1. Thus we have S3/2/C — — 3, which is totally (anti)-correlated to the curvature 
perturbation and already contradicts with the current observation [1]. This shows that 
a scenario with thermally produced gravitino DM cannot be viable in the modulated 
reheating scenario, in which large non-Gaussianity can be generated, because of too large 
isocurvature fluctuations. 

In fact, a similar argument also holds when gravitinos are produced non-thermally 
from the decay of some heavy scalar field such as the inflaton or moduli in the case they 
acquire a nonvanishing vacuum expectation value [15], which dominates over the other 
channels [18]. In this case, the yield can be written as 

^ = l B ^> (27 > 

s I m 3 /2 

where -B 3 / 2 is the branching ratio of the decay into gravitinos. When gravitinos are pro- 
duced from the jets, B 3 / 2 should be understood as those including its multiplicity. Since 
the reheating temperature is related to the decay rate of the inflaton as Tr oc T 1//2 , then 
£>3/ 2 oc T^ 2 , we have the T^-dependence of the n 3 / 2 /s, in the case of non-thermal produc- 
tion, as 

=5* cc * (28) 
s T R 

Thus isocurvature fluctuations are also generated when gravitinos are produced from some 
heavy scalar decay, which can be written as 

S3/2 = (29) 
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which is different in sign but the same in size as in the case that gravitinos are produced 
thermally. Note that this isocurvature fluctuation has totally positive correlation to the 
curvature perturbation. Hence if DM consists of gravitinos, large isocurvature fluctuations 
are generically produced in the modulated reheating scenario, which is inconsistent with 
cosmological observations such as CMB. 

There is another well-known mechanism to generate large local-type non-Gaussianity, 
called the curvaton mechanism. In this scenario, gravitinos are produced from the decays 
of both the inflaton and the curvaton. As discussed in Ref. [11, 19], in the former case, 
gravitino DM is disfavored because of too large isocurvature fluctuations. On the other 
hand, in the latter case, gravitino DM is permitted only when the curvaton almost domi- 
nates the energy density of the universe to avoid large isocurvature fluctuations, in which 
the non-Gaussianity is small. Then, it is concluded that gravitino DM is disfavored also 



in the curvaton scenario if the large non-Gaussianity is observed [13 Therefore, gravitino 
dark matter is disfavored once large local-type non-Gaussianity of primordial fluctuations 
is observed. 

Some comments are in order. First gravitinos can also be produced by the decay of the 
next-to-the-lightest supersymmetric particle (NLSP). However, the detailed calculations 
show that the constraints from BBN are severer than those from the overclosure of the 
universe, irrespective of the particle kind of NLSP as long as it is the MSSM particle [21]. 
Therefore, gravitino DM can be realized only for the case that most of gravitinos should 
be produced much earlier, either by the scattering processes during reheating or by the 
non-thermal decay of heavy scalars, both of which we paid attention to in this Letter. 

Second, a similar argument also holds for the case where LSPs, which are DM, are pro- 
duced from the decay of gravitinos. In this case, LSPs carry large isocurvature fluctuations 
originating from gravitinos, which again contradicts with the present constraint. Thus, 
DMs originating from gravitinos are also disfavored once large local-type non-Gaussianity 
is detected^!. 

Finally we comment on the case of axino DM. Axinos are produced by thermal scat- 
tering during reheating and their abundance also depends on the reheating tempera- 
ture [22-25]. Therefore, as with the case of gravitinos, if such axinos become DM, they also 
acquire large isocurvature fluctuations in the modulated reheating scenario, which again 
contradicts with the present observational constraints. Unlike the gravitino dark matter, 
however, non-thermally produced axinos by the decay of NLSP [24-26] or Q-balls [27] can 
be dominant with avoiding the BBN constraints. For such cases, isocurvature constraints 
can be waived. 



Acknowledgments: 

We thank Takahiro Tanaka for the collaboration at the early stage. We are also grateful 
to Takashi Hamazaki, Masahiro Kawasaki, Kazunori Kohri, and Fuminobu Takahashi for 

# 4 For the case of baryonic isocurvature fluctuations, sec Ref. [20]. 

# 5 We thank M. Kawasaki and F. Takahashi for pointing out this possibility. 



7 



useful discussions. This work is supported by JSPS Grant-in-Aid for Scientific research, 
No. 19740145 (T.T.), No. 19340054 (J.Y.), and No. 21740187 (M.Y.). S.Y. is supported in 
part by Grant-in-Aid for Scientific Research on Priority Areas No. 467 "Probing the Dark 
Energy through an Extremely Wide and Deep Survey with Subaru Telescope". He also 
acknowledges the support from the Grand-in-Aid for the Global COE Program "Quest for 
Fundamental Principles in the Universe: from Particles to the Solar System and the Cos- 
mos " from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) 
of Japan. We would like to thank the organizers of the IPMU workshop on "Focus week on 
non-Gaussianities in the sky" and the GCOE/YITP workshop YITP-W-09-01 on "Non- 
linear cosmological perturbations" for their hospitality, during which a part of this work 
is done. 



References 

[1] E. Komatsu et al. [WMAP Collaboration] , [astro-ph]. 

[2] K. M. Smith, L. Senatore and M. Zaldarriaga. !arXiv:0901.2572l [astro-ph], 

[3] A.H. Guth, Phys. Rev. D23, 347 (1981); K. Sato, Mon. Not. Roy. Astron. Soc. 
195, 467 (1981); A. A. Starobinsky Phys. Lett. 91B, 99 (1980); A. D. Linde, Phys. 
Lett. 108B, 389 (1982); A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, 1220 
(1982); A.D. Linde, Phys. Lett. B 129, 177 (1983). 

[4] S.W. Hawking, Phys. Lett. 115B, 295 (1982); A.A. Starobinsky, Phys. Lett. 117B, 
175 (1982); A.H. Guth and S-Y. Pi, Phys. Rev. Lett. 49, 1110 (1982). 

[5] G. Dvali, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69, 023505 (2004) 
|arXiv:astro-ph/030359l] . 

L. Kofman, arXi v: astro-ph /0303614[ 

[6] S. Mollerach, Phys. Rev. D 42, 313 (1990). 

A. D. Linde and V . F. Mukhanov, Phys. Rev. D 56, 535 (1997) 
|arXiv:astro-ph/9610219| . 

[7] K. Enqvist and M. S. Sloth, Nucl. Phys. B 626, 395 (2002) [arXiv:hep-ph/0109214| ; 

D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) |arXiv:hep-ph/0110002| ; 

T. Moroi and T. Takahashi, Ph ys. Lett. B 522, 215 (2001) [Erratum-ibid. B 539, 303 
(2002)] |arXiv:hep-ph/01 10096] . 

[8] E. W. Kolb, A. Riott o and A. Vallinotto, Phys. Rev. D 73, 023522 (2006) 
|arXiv:astro-ph/05liT98| . 

D. H. Lyth, JCAP 0511, 006 (2005) [arXiv: astro-ph/05 10443| . 



8 



L. Alabidi and D. Lyth, JCAP 0608, 006 (2006) [arXiv:astro-ph/0604569] . 
M. Sasaki, Prog. Theor. Phys. 120, 159 (2008) |arXiv:0805.0"974 [astro-ph]]. 

[9] M. Zaldarriaga, Phys. Rev. D 69, 043508 (2004) [arXiv: astro-ph/0306006| . 

A. Mazumdar, Phys. Rev. Lett. 92, 241301 (2004) |arXiv:hep-ph/0306026| . 

K. Ichikawa, T. Suyama, T. Takahashi and M. Yamaguchi, Phys. Rev. D 78, 063545 
(2008) [arXiv:0807.3988l [astro-ph]]. 

[10] T. Suyama and M. Yamaguchi, Phys. Rev. D 77, 023505 (2008) [ arXiv:0709.2545l 
[astro-ph]]. 

[11] D. H. Lyth, C. Ungar elli and D. Wands, Phys. Rev. D 67, 023503 (2003) 
|arXiv:astro-ph/0208~0~55l . 

D. H. Lyth and D. Wands, Phys. Rev. D 68, 103516 (2003) |arXiv:astro-ph/0306500] . 

[12] N. Bartolo, S. Matar rese and A. Riotto, Phys. Rev. D 69, 043503 (2004) 
[arXiv:hep-ph/0309033| . 

K. Enqvist and S. Nurmi, JCAP 0510, 013 (2005) |arXiv:astro-ph/0508573] . 

K. A. Malik and D. H. Lyth, JCAP 0609, 008 (2006) [arXiv: astro-ph/060438 7| . 

M. Sasaki, J. Valivii ta and D. Wands, Phys. Rev. D 74, 103003 (2006) 
|arXiv:astro-ph/0607627| . 

Q. G. Huang. la7XTv:0801.0467l [hep-th]. 

K. Ichikawa, T. Suyama, T. Takahashi and M. Yamaguchi, Phys. Rev. D 78, 023513 
(2008) [arXiv: 0802 .4138 [astro-ph]]. 

K. Enqvist and T. Takahashi, JCAP 0809, 012 (2008) [arXiv:0807.3069l [astro-ph]]. 

Q. G. Huang, JCAP 0811, 005 (2008) (arXiv: 0808 .17931 [hep-th]]. 

Q. G. Huang and Y. Wang, JCAP 0809, 025 (2008) [arXiv: 0808 .11681 [hep-th]]. 

[13] M. Alishahiha, E. Sil verstein and D. Tong, Phys. Rev. D 70, 123505 (2004) 
|arXiv:hep-th/0404084| . 

N. Arkani-Hamed, P. Creminelli, S. Mukohyama and M. Zaldarriaga, JCAP 0404, 
001 (2004) |arXiv:hep-th/0312100| . 

[14] J. R. Ellis, J. E. Kim and D. V. Nanopoulos, Phys. Lett. B 145, 181 (1984). 

T. Moroi, H. Murayama and M. Yamaguchi, Phys. Lett. B 303, 289 (1993). 

M. Bolz, A. Brandenburg and W. Buchmuller, Nucl. P hys. B 606, 518 (2001) 
[Erratum-ibid. B 790, 336 (2008)] [arXiv:hep-ph/0012052] . 

J. Pradler and F. D. Steffen, Phys. Rev. D 75, 023509 (2007) |arXiv:hep-ph/0608344] . 



9 



[15] S. Nakamura and M. Yamaguchi, Phys. Lett. B 638, 389 (2006) 
|arXiv:hep-ph/0602d8i| . 

T. Asaka, S. Nakamu ra and M. Yamaguchi, Phys. Rev. D 74, 023520 (2006) 
|arXiv:hep-ph/0604132] . 

M. Kawasaki, F. Takahashi and T. T. Yanagida, Phys. Lett. B 638, 



8 (2006)|arXiv:hep-ph/0603265'|; Phys. Rev. D 74, 043519 (2006) 
|arXiv:hep-ph/0605297] ~ 

M. Endo, K. Hamagu chi and F. Takahashi, Phys. Rev. Lett. 96, 211301 (2006) 
|arXiv:hep-ph/0602d6i| . 

M. Endo, F. Takahashi and T. T. Yanagida, Phys. Rev. D 76, 083509 (2007) 
|arXiv:0706.0986l [hep-ph]]. 

[16] C. Hikage, K. Koyama, T. Matsubara, T. Takahashi and M. Yamaguchi, 
larXiv:0812.3500l [astro-ph]. 

[17] A. A. Starobinsky, JETP Lett. 42 (1985) 152 [Pisma Zh. Eksp. Teor. Fiz. 42 (1985) 
124]. 

Y. Nambu and A. Taruya, Class. Quant. Grav. 13, 705 (1996) 
|arXiv:astro-ph/94li0T3| . 

M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 (1996). 
M. Sasaki and T. Tanaka, Prog. Theor. Phys. 99, 763 (1998). 

H. Kodama and T. Hamazaki, Phys. Rev. D 57, 7177 (1998) |arXiv:gr-qc/9712045] . 
D. H. Lyth, K. A. Malik, and M. Sasaki, J. Cosmol. Astropart. Phys. 05, 004 (2005). 
T. Hamazaki, Phys. Rev. D 78, 103513 (2008) [arXiv:0811.2"366l [astro-ph]]. 

[18] M. Hashimoto, K. I. Izawa, M. Yamaguchi and T. Yanagida, Prog. Theor. Phys. 100, 



395 (1998) |arXiv:hep-p h/9804411|; 

K. Kohri, M. Yamagu chi and J. Yokoyama, Phys. Rev. D 70, 043522 (2004) 
|arXiv:hep-ph/0403043| . 

[19] M. Lemoine, J. Martin and J. Yokoyama, larXiv:0903.5428l [astro-ph. CO]; 
larXiv:0904.0T26l [astro-ph. CO]. 

[20] T. Moroi and T. Takahashi, Phys. Lett. B 671, 339 (2009) |arXiv:0810.0T89l [hep-ph]]. 

[21] M. Kawasaki, K. Kohri, T. Moroi and A. Yotsuyanagi, Phys. Rev. D 78, 065011 
(2008) [arXiv:0804.3745l [hep-ph]]. 

[22] K. Rajagopal, M. S. Turner and F. Wilczek, Nucl. Phys. B 358, 447 (1991). 



10 



[23] L. Covi, L. Roszk owski and M. Small, JHEP 0207, 023 (2002) 
|arXiv:hep-ph/0206119| . 

[24] L. Covi, H. B. Kim, J. E. Kim and L. Roszkowski, JHEP 0105, 033 (2001) 
|arXiv:hep-ph/0101009| . 

[25] A. Brandenburg and F. D. Steffen, JCAP 0408, 008 (2004) [arXiv:hep-ph/0405158] . 
A. Freitas, F. D. Steffen, N. Tajuddin and D. Wvler. larXTv:0904.3218l [hep-ph]. 

[26] L. Covi, J. E. Kim and L. Roszkowski, Phys. Rev. Lett. 82, 4180 (1999) 
|arXiv:hep-ph/9905212] . 

[27] L. Roszkowski and O. Seto, Phys. Rev. Lett. 98, 161304 (2007) 
|arXiv:hep-ph/06080"T3] . 

O. Seto and M. Yamaguchi, Phys. Rev. D 75, 123506 (2007) |arXiv:0704.05l0l [hep- 
ph]]. 



11 



